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Abstract : The pion structure function is studied in the Nambu and Jona-Lasinio (NJL) 
model. We calculate the forward scattering amplitude of a virtual photon from a pion 
target in the Bjorken limit, and obtain valence quark distributions of the pion at the low 
energy hadronic scale, where the NJL model is supposed to work. The calculated distribu- 
tion functions are evolved to the experimental momentum scale using the Altarelli-Parisi 
equation. The resulting distributions are in a reasonable agreement with experiment. We 
calculate also the kaon structure function and compare the ratio of kaon to pion valence 
u-quark distributions with experiment. 
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Recent experimental observations of deep inelastic scattering (DIS) provide us with a 
large amount of data about the structure of hadrons. These results |], || indicate that our 
understanding of hadron structure is still incomplete. Quantum chromodynamics (QCD) 
is believed to be the theory of the strong interaction, and perturbative QCD is consistent 
to the measured Q 2 dependence of the structure function However, QCD is not able to 
predict directly the structure function itself, since non-perturbative effects of confinement 
and dynamical symmetry breaking are not yet understood in the context of QCD. On 
the other hand, phenomenological quark models, which possess some non-perturbative 
features of QCD, give remarkable predictions of static hadron properties such as masses 
and magnetic moments. Hence, it is of our great interest to clarify a connection between 
the DIS phenomena and the low energy non-perturbative features in terms of such models 
as effective theories of QCD. Such a connection would enable us to use the DIS data to 
constrain the models of low energy QCD. 

Recently, the nucleon structure functions are discussed in terms of the MIT bag mode 
0, U as well as other models ||. Those results are based on the assumption that the 
structure functions at low energy hadronic scale Q 2 = Qq (unknown) are obtained by 
calculating the twist 2 matrix elements within the effective models, since the twist 2 
operators are dominant in the Bjorken limit and higher twist terms vanish at high enough 
Q 2 [|J. Once one obtains them at the hadronic scale, these parton distributions are 
evolved to the experimental scale by using the perturbative QCD, and the comparison 
with experiment can be made. 

In this letter, we concentrate on the pion structure function. The pion is believed to be 
the Goldstone boson due to the spontaneous breakdown of the chiral symmetry (SB^S). 
Thus, the information of the pion structure function gives us deeper understandings of 
SB^S. The experimental data of the pion structure function are extracted from 7rN scat- 
tering using the Drell-Yan process pj. The valence, sea quark, and gluon distribution 
functions are obtained from a recent re-analysis of these data ||. 

We use the Nambu and Jona-Lasinio (NJL) model ||, where the chiral invariance is 
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the main ingredient. The NJL model describes the pion as a collective qq excitation of the 
non-perturbative QCD vacuum Recently, the generalized SU(3) / NJL model with the 
U(l)yi anomalous term is demonstrated to reproduce the meson properties successfully, 



in spite of the lack of confinement |10|, fy], |12| . This model is also applied to the chiral 
phase transition at finite temperature and density [13]. All these results indicate that the 
NJL model possesses the essential features of QCD dynamics. 



The SU(3) / NJL lagrangian is given by [|T( 



0,11 



£njl = ip{iYd^ - m)ip 

+G s [{i,\ l i,) 2 + (^7sA^) 2 ] - GV[(?HAV0 2 + (^7sA^) 2 ] , ( ■'■ ) 



where ip denotes the quark field with current mass m, \ are SU(3) flavor matrices, and 
Gs, Gy are the coupling constants. In this model, the quarks acquire the constituent 
masses dynamically due to SB%S. The meson masses are obtained by solving the Bethe- 
Salpeter equation. This model describes the SU(3)/ meson properties very well with the 



parameters fixed by the pion and kaon properties |T0| , |TT| , [12 1 • Using this lagrangian, we 
will calculate the parton distribution of pion with no free parameter. 

Unfortunately, the NJL model is a non-renormalizable theory, and this model requires 
the finite momentum cutoff A ~ 1 GeV, which is identified with a scale of SB%S. One may 
understand the physical meaning of the cutoff as an approximate realization of "asymp- 
totic freedom" in the NJL model. This means that the interaction between two quarks 
with the relative momentum larger than A is turned off, and two particles are free in 
such a high momentum scale. Hence, the structure functions obtained in the NJL model 



exhibit the Bjorken scaling in the deep inelastic limit [JTTJ] . If we use the cutoff in a Lorentz 
invariant manner, we can calculate the quark distributions as a function of the Bjorken 
x. 

As mentioned above, we would like to calculate the structure function in the NJL 
model at a boundary to the QCD perturbation. The hadronic tensor is written by 



the structure functions F\ and F 2 in the scaling limit; 



+ —(P^-^^)(P^-^9u)F 2 {x), (2) 



where 



F a (z) = z£eft ft (s) +*(*)], F 1 (x) = ^. 

and are the momentum distributions of the i-flavor quark and antiquark. The 
hadronic tensor is related to the forward scattering amplitude through the optical 
theorem ||. 



= ^ImT^ . (3) 

Thus, we evaluate T^ v in the NJL model to get the structure functions. The forward 
scattering amplitude of a virtual photon with q from a pion is defined by 



T ia , = iJ *te*<p\T[Jp®J v (p)]\p)o , (4) 

where \p) is a pion state of momentum p^ , and the momentum delivered by the current. 
The current is J M = ijj'y^Qijj with Q being the charge operator. We compute the "handbag 
diagrams" for T^ Vl which is illustrated in fig. la and lb. The matrix element (flj) in the 
NJL model is given by 



d 4 k 1 

Tp, = i I j^^Tt^^Q^uQT-] + T+ term , (5) 



where 



T_ = S F (k- q)g 7Tqq r + i'y 5 S F (k -q- p)g 7Tqq r^i'y 5 S F (k - q) . (6) 

Here, g nqq is the coupling constant among two quarks and a pion obtained in the NJL 
model |ll|, |I2|| , and r± = (ti ± ir 2 )/\/2- T- and T + denote scattering of a pion from 
quark and antiquark, respectively, where T + has a similar expression as (|6]). Note that 
the quark propagator contained in T_ has the dynamically generated mass obtained by 
solving the gap equation || [L(|]. Thus, S F (k) is written as 



S F (k) ' 



where M is the constituent quark mass. 

The structure function obtained by the calculation of the "handbag diagram", shown 
in figs. la and lb, coincides with the leading twist contribution of the operator product 



expansion in the Bjorken limit [[14]]. Therefore, these diagrams are enough for our purpose, 
because the structure functions are dominated by the twist 2 contributions in the DIS 
limit. Higher twist terms may be negligible at sufficiently high Q 2 H. 

Following the work of Landshoff et al. [j[5f| , we will carry out the integration of (|^) in 
the Bjorken limit; 



n2 2 v-q Q 2 n A 

Q = — q — > oo , v = > oo , x = : faxed . 

Here x is the so-called Bjorken x, and the pion mass. We introduce the integral 
variables 



K = zp^ + yq^ + k m , (7) 
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where satisfies k ■ p = k • q = 0. Thus, is spacelike (k 2 < 0) and is effectively two 
dimensional. Calculating the traces and neglecting irrelevant terms in the Bjorken limit 
15 1, we obtain 



Sir 2z 1 
T ^ = nTT^i Nc &m / dzdyd 2 K(-g^ v + p^p v + (p M g„ + 



9(2vr) 

x[t 1 ( / u 2 ,s) + zt 2 (/i 2 ,s)] 1 



2 — x ' 



where 



ti(/i 2 ,s) 
t 2 (/x 2 ,s) 



1 



V-M 2 ) 



(/i 2 — M 2 + is) 2 s- M 2 + ie 

( 8 — M 2 — v 2 ) 

(jj 2 — M 2 + is) 2 s — M 2 + is P} 



and 



(k — q) 2 = zy + 



s = (k — q — p) 2 = (z — l)(y — m 2 ) + k 2 



We change the variable y = Im^viy — 1) + zp |]15 |, and // and s are the invariant masses 



of the struck quark and spectator. Performing the z- integral in eq.flq), we find [15, I I 



nn 8 1 »r 2 r 2 ^ 1 / 

x J dyd 2 ; /t[ti(/i 2 , s) + xt 2 (/i 2 , s)] . (9) 

Consider now the y-integral in the complex y-plane. The s-propagator has a pole at 
y = (k 2 — M 2 )/(l — x) + m 2 + is /{I — x), and the //-propagator has a double pole at 



y = (M 2 — k 2 )/x — ie/x. For x > 1 or x < 0, all these singularities occur on the same 
side of the real y-axis, and @ gives a zero result |L5], However, if < x < 1, the 
integration over y no longer vanishes. Integrating @ by y and taking the imaginary 
part, we can get the quark distribution by use of the optical theorem (|3]). We change the 
integral variable k to fi 2 , we find the following expression for the quark distribution: 



x#(m 2 x(l -x) -xM 2 - (1 -x)fi 2 ) . (10) 



Here, 6 is the usual step function which arises from the spacelike condition of k. It is easily 
shown that the structure functions obtained from (pi]) exhibit the Bjorken scaling JTTJ] . We 



identify flTU| ) as the valence quark distribution of the pion qv{x), with the normalization 
J dxqv(x) = 1. The integral (0) goes to infinity due to the non-renormalizability of the 
NJL model. Thus, we shall use the momentum cutoff. There exists an ambiguity in the 
introduction of the NJL cutoff for DIS phenomena |16|, [XT]. We introduce the Euclidean 
momentum, t 2 E = —fi 2 + m 2 x — xM 2 /(l — x), and simply assume that the four momentum 
of the struck quark has the sharp Euclidean cutoff, #(A 2 — t 2 E ), which is consistent to the 
NJL cutoff. We examine also the exponential cutoff used in ref. |TB[ and other plausible 



forms for completeness. We obtain similar results with ours except for the behavior of 
q(x) around x ~ 1 which depends on the cutoff scheme. 

Integrating ( PH| ) by £e(jU 2 ) with the cutoff, we obtain the valence distribution of pion at 
the low energy hadronic scale, which is shown in fig.2. Here, we use the NJL parameters 
fixed by the pseudoscalar meson properties A = 900MeV, G5A 2 = 2GyA 2 = 4.67, m u = 



vrid = 5.5MeV and m s = 135MeV |L0], |TT, With this parameter set, we obtain the 
constituent masses, M u = = 350MeV and M s = 53AM eV. Our results on the quark 
distribution depend only slightly on the choice of the parameters within the reasonable 
parameter ranges, once these parameters are constrained by the pion mass and the decay 
constant. Note that the resulting distribution shows a correct behavior around x ~ 1 



; qv{ x ) ~~ > as x — > 1 |TT|. It is the advantage of our use of the NJL model, which is 
Lorentz invariant, as compared with the MIT bag model where the translational invariance 
is broken ||. We also note that this low energy scale structure function has no physical 
meaning at this scale, since "real" structure function at the low energy scale receives 
non- negligible contributions from all twist operators. Our calculated result in fig. 2 plays 
a role of only the boundary condition of the structure function for the Q 2 evolution. 
We take the low energy hadronic scale at Qq = (0.5Ge\^) 2 , which is used in ref. [|18 



At this scale, the running coupling constant is still small; a s (Ql)/ii ~ 0.3. Indeed, the 
inclusion of the second order QCD corrections gives a small change for the Q 2 evolution 
from our result within 10% []18|| . We may understand intuitively the physical implication 
of this scale Ql as compared to the valence quark core radius of the pion (r 2 ) core , noted 
by Brown et al. |L9]]. Their value is consistent to our low energy scale: 



(r 2 ) co , e ~ (0.35fm) 2 ~ l/{0.5GeVf 



We use the first order Altarelli-Parisi equation [3| for the Q 2 evolution of valence 
distributions with Aqcd = 250MeV. We show in fig. 2 (solid curve) the result at Q 2 = 
20GeV 2 , which is compared with the experimental curve (dash-dotted curve) extracted 
from the Drell-Yan process ||. We find a reasonable agreement with experiment. If 
we vary the low energy scale Ql to a larger value Ql = (0.75GeV) 2 , the change of 
our distribution is of order 20distribution function at the experimental scale, it seems 
necessary to choose Ql as small as Ql = (0.5Ge^) 2 . We compare the pion structure 
function F n (x) with the experiment in fig. 3. It shows a good agreements for x > 0.3. In 
the low x region, the structure function is dominated by the sea quark distribution ||. 
We ought to calculate the sea quark contributions |L6|, [L7j with the next-to-leading order 
QCD perturbation, which is under consideration. 

We also calculate the first two moments of the valence quark distributions, (x) and 



s 



(x 2 ), which are defined as 



(x) = / dxxq v (x) , (x 2 ) = / dxx 2 q v (x) . 
Jo Jo 



The resulting values are shown in table 1 with the experimental values M and the results 



of the lattice QCD calculations Our results reproduce the data very well. 

Finally, we discuss the kaon structure function. We can extract the effects of the SU(3) 
flavor symmetry breaking from quark distributions of the kaon. The valence quark dis- 
tributions are obtained in a similar manner using the NJL model. The u- and s-quark 
distributions of the kaon at the low momentum scale and those at Q 2 = 20GeV 2 are shown 
in fig.4. Our results indicate that the heavy s-quark carries a larger fraction of kaon mo- 
mentum than the light u(d) quark, as expected. The resulting distributions show a similar 



behavior with the previous results obtained by the analysis of the Regge trajectories |22 



We also show in fig. 5 the ratio of kaon to pion valence u-quark distributions uk/u^ at 
Q 2 = 20GeV 2 with experimental data [p3fl . This ratio is sensitive to the mass difference 
of the constituent quark masses. In fact, our calculations give a simple expression of this 
ratio around x ~ 1 at the low momentum scale as 



~ (M u /M s ) 2 ~ 0.5 . 

This value is quite reasonable in view of the experimental findings. 

In conclusion, we have studied the pion structure function at the DIS scale using the 
NJL model as a low energy effective theory of QCD. The resulting distribution at the low 
momentum scale is evolved to the experimental scale, at which it compares with the DIS 
data very well. We have also studied the kaon structure function. Reflecting the SU(3) 
flavor symmetry breaking due to the bare current mass difference, the x dependence of the 
ratio of kaon to pion valence u-quark distributions is found consistent with experiment. 
Our results indicate that the low energy effective quark model is able to describe the DIS 
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phenomena. We shall point out that the sea quark degrees of freedom can be naturally 
included in the framework of the NJL model, which is needed to compare with the pion 
structure function directly. It is possible to calculate the sea quark distributions in the 
NJL model as higher loop corrections |16|, |ITf . 

We are grateful to T. Hatsuda and T. Kobayashi for illuminating discussions on the 
deep inelastic phenomena. 
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Table Caption 



Table 1 : The first two moments of the pion distribution at Q 2 = 49GeV 2 . The theoret- 
ical calculations of the NJL model are shown in the second row, and compared with 
the experimental values || (third row) and the lattice QCD results |2]J (last row). 



Table 1 





2(x) 


2(x 2 ) 


Theory 


0.41 


0.16 


Experiment 


0.40±0.02 


0.16±0.01 


Lattice 


0.46±0.07 


0.18±0.05 
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Figure Captions 



Fig. 1 : The forward Compton scattering amplitude ("handbag diagram") of the pion. 
The solid line represents the quark. The pion and the virtual photon are depicted 
by the dashed and wavy lines. For details and notation, see text. 

Fig. 2 : The valence quark distributions of the pion at the low energy scale Q 2 = Qq 
(dotted curve) and at Q 2 = 20GeV 2 (solid curve) as a function of the Bjorken x. 
The experimental fit |§ is depicted by the dash-dotted curve. 

Fig. 3 : The pion structure function at Q 2 = 20GeV 2 . The experimental data are taken 
from the NA3 experiment |7J . The theoretical prediction of the N JL model is depicted 
by the solid curve. Here, we use F n (x) = Kxqv(x) with the sea quark distributions 
being neglected, where we take the K- factor, K—1.5 ||. 

Fig. 4 : The valence quark distributions of the kaon at low energy scale and at Q 2 = 
20GeV 2 in the NJL model. The u- and s-quark distribution at Qq are depicted by 
the dotted and dashed curves. The solid and dashed-dotted curves represent the u- 
and s-quark valence distributions at Q 2 = 20GeV 2 . 

Fig. 5 : The ratio of kaon to pion valence u-quark distributions 

UK/u n at Q 2 = 20GeV 2 scale. The theoretical result is depicted by the solid curve. 
The experimental values with error bars are taken from ref. 
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Fig. 1 
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